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Abstract 

We extend the analysis of Higgs-regulated planar amplitudes of A/" = 4 supersymmet- 
ric Yang-Mills theory to four loops for the four-gluon amplitude and to two loops for the 
five-gluon amplitude. Our calculations are consistent with a proposed all-loop ansatz for 
planar MHV n-gluon amplitudes that is the analog of the BDS ansatz in dimensional regu- 
larization. In all cases considered, we have verified that the IR-finite parts of the logarithm 
of the amplitudes have the same dependence on kinematic variables as the corresponding 
functions in dimensionally-regulated amplitudes (up to overall additive constants, which we 
determine) . 

We also study various Regge limits of A/" = 4 SYM planar n-gluon amplitudes. Euclidean 
Regge limits of Higgs-regulated n > 4 amplitudes yield results similar in form to those found 
using dimensional regularization, but with different expressions for the gluon trajectory 
and Regge vertices resulting from the different regulator scheme. We also show that the 
Regge limit of the four-gluon amplitude is dominated at next-to-leading-log order by vertical 
ladder diagrams together with the class of vertical ladder diagrams with a single H-shaped 
insertion. 
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1 Introduction 



The assumption of dual conformal symmetry has proven useful for understanding the structure 
of the large- limit of higher-loop scattering amplitudes in TV = 4 supersymmetric Yang-Mills 
theory. Dual conformal invariance characterizes [1] the set of scalar diagrams that forms a basis 
for the computation of higher-loop amplitudes in the unitarity-based approach of refs. pHH]. 
Moreover, coupled with the assumption of Wilson loop/MHV scattering amplitude duahty [7HTU]. 
dual conformal symmetry has been used fXl] to explain why maximally-helicity-violating planar 
n-gluon amplitudes obey the Bern-Dixon-Smirnov (BDS) ansatz [12] for n = 4 and 5. It also 
helps to explain why the BDS ansatz fails to predict the scattering amplitude starting at n = 6, 
as observed in refs. [T3HTB] . and constrains the form of the discrepancy, called the remainder 
function (for reviews see refs. [IZIITB]). Wilson loop/MHV amphtude duality has been used to 
determine the precise form of this remainder function for n = 6 in refs. p ^[T9]42T| . 

Scattering amplitudes in massless gauge theories possess infrared (IR) divergences that must 
be regulated. It is desirable that the regulator preserve as many symmetries as possible. While 
dimensional regularization explicitly breaks dual conformal symmetry (see, e.g., ref. an 
alternative Higgs regulator proposed by Alday, Henn, Plefka, and Schuster [23] leaves the dual 
conformal symmetry unbroken. In this approach, the supersymmetric Yang-Mills (SYM) theory 
is considered on the Coulomb branch where scalar vevs break the gauge symmetry, causing some 
of the gauge bosons to become massive through the Higgs mechanism. Planar gluon scattering 
amplitudes on this branch can be computed using scalar diagrams in which some of the internal 
and external states are massive, regulating the IR divergences of the scattering amplitudes. (For 
earlier applications of a massive IR regulator, see refs. [7t l2^H27] : low-energy amplitudes in the 
Higgsed phase of SYM have been studied in ref. [28].) The diagrams remain dual conformal 
invariant, however, provided that the dual conformal generators are taken to act on the masses 
as well as on the kinematical variables, a generalization referred to as extended dual conformal 
symmetry. 

The assumption of extended dual conformal symmetrjfl severely restricts the number of di- 
agrams that can appear in the scattering amplitudes. In particular it forbids loop integrals 
containing triangles, which have indeed recently been shown to be absent from one-loop ampli- 
tudes on the Coulomb branch [29] (a related discussion was given in ref. [26]). At one point in 
moduli space, all the lines along the periphery of the diagrams have mass m, while the external 
states and the lines in the interior of the diagram are all massless. This is believed to be sufficient 
to regulate all IR divergences of planar scattering amplitudes!! The original SYM theory is then 
recovered by taking m small. 

Using the Higgs regulator described above, the J\f = 4 SYM planar four-gluon amplitude was 
computed at one and two loops in ref. ^3] and at three loops in ref. [50], assuming that only 
integrals invariant under extended dual conformal symmetry contribute. It was shown that, at 

^or rather the slightly stronger assumption (which holds in every known example) that every amplitude can 
be expressed in an integral basis in which each element is dual conformal invariant 

^More precisely, diagrams which cannot be rendered finite in this manner also cannot be rendered finite in 
off-shell regularization, and it is believed that such diagrams never contribute to the amplitude [HI. 
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least through three loops, the Higgs-regulated four-gluon amplitude obeys an exponential ansatz 
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analogous to the BDS ansatz in dimensional regularization. Here M4^{s,t) is the ratio of the all- 
orders planar amplitude to the tree-level amplitude, 7(a) is the cusp anomalous dimension [31], 
and QoiO') and 64(0) are analogs of functions appearing in the BDS ansatz. One of the advantages 
of the Higgs-regulated ansatz (11. ip is that IR divergences take the form of logarithms of m^; con- 
sequently, the L-loop amplitude may be computed by simply exponentiating logM4(s,t) without 
regard for the 0{m^) terms since they continue to vanish as m — > even when multiplied by log- 
arithms of m?. Putting it another way, in order to test eq. (11 .11) one need not compute any 0{m'^) 
terms of the Higgs-regulated L-loop amplitudes because they cannot make any contribution to 
the IR-finite part of logM(s,t). This stands in stark contrast to dimensional regularization, 
where checking the BDS ansatz at each additional loop order requires recalculation of two more 
terms in the e expansion of each lower loop integral. 

One of the results of the current paper is to extend the computation of the Higgs-regulated 
four-gluon amplitude to four loops. Following the observed behavior through three loops, we 
assume that the Higgs-regulated four-loop amplitude can be expressed as the same linear com- 
bination of eight scalar integrals as the dimensionally-regulated amplitude. Specializing to the 
kinematic point s = t, we demonstrate that the result is consistent with the exponentiation 
of IR divergences built into the BDS ansatz (11. ip . and we confirm, with significantly improved 
numerical precision, the value of the four- loop cusp anomalous dimension found in refs. [321E3] • 

Taking the Regge limit s 3> t, the ansatz (II. ID implies that the Higgs-regulated four-gluon 
amphtude exhibits exact Regge behavior [Sj lSlHSU] 



where the all- loop-orders Regge trajectory is 
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The coefficient of log t, the cusp anomalous dimension, is independent of the IR regulator, while 
the constant part is scheme-dependent. In ref. ^D], we verified that the Higgs-regulated four- 
gluon amplitude obeys eq. (II. 2p through three-loop order, and in the current paper, we extend 
this to four loops. 

Equation (II. ip implies that the L-loop amplitude has leading log (LL) expansion 



M, 



:-ii 



log' 



+ 



+ 



L! 

(-1) 



t 



log' 



TT 



TT 



2(L 

s 



6(L-2)! 



log 



L-l 



(L-2) 



log 



L-2 



log 



L-l 



;i-4) 



3 



where the LL log'^ s term depends only on the lowest-order term 7(a) = —4a + 0{a'^) of the cusp 
anomalous dimension, while the next-to-leading log (NLL) log'^~^ s term depends on the 0{a'^) 
terms of 7(a) and Qo{0')- We showed in ref. [30] that the LL term stems entirely from a single 
scalar diagram, the vertical ladder, in the Regge limit of the Higgs-regulated loop expansion]^ In 
this paper, we show that the NLL contribution to the L-loop amplitude is given by the subleading 
term of the vertical ladder, together with the leading contributions from a set of (L — 1) diagrams, 
consisting of vertical ladder diagrams with a single H-shaped insertion. We explicitly confirm 
this through five loops via a Mellin-Barnes computation, and present an argument (subject to 
certain reasonable assumptions) for its validity to all loops. 

A second thrust of this paper is to consider higher-point amplitudes in Higgs regularization. 
We compute the Higgs-regulated five-gluon amplitude at one- and two-loops and establish the 
iterative relation 
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is the IR-finite part of the one-loop amplitude, which is the same as in dimensional regularization, 
up to an additive constant. We argue that the parity-odd part of M^^^ is at most O(m^) in Higgs 
regularization. 

Based on the iterative relation (11.51) . we propose that the generalization of eq. (II. ip to the 
planar MHV n-gluon amplitude in Higgs regularization takes the form 
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and TZn vanishes for n = 4 and n = 5. For n > 6, we expect the remainder function TZn to be 
equal to its counterpart in dimensional regularization. 

Various Regge limits of the n-gluon amplitude forn > 5 can be defined (see, e.g., refs. [iOlHT] ). 
We consider the single and double Regge limits (in the Euclidean region) of the Higgs-regulated 
five-gluon amplitude up to two loops. In both of these cases, the double logarithm in logMs 
cancels out, leaving single logarithmic dependence on the large kinematic variable. 



^ There are two ways of taking the Regge hmit of a Higgs-regulated amphtude. One can either (a) first take 
the limit ^ s,t, and then s ^ i, or (b) first take the limit s ^ t, rri^, and then <C t. We demonstrated 
in ref. [30] that the amplitude is independent of the order of limits, at least through three loops. The Regge 
behavior of individual diagrams, however, can depend on the order in which the limits are taken. The dominance 
of vertical ladder diagrams is only valid in the Regge (b) limit, and the discussion in the text assumes this order 
of limits. 
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In ref. pJD], an alternative approach to the Regge hmit for four-gluon amphtudes was con- 
sidered by taking a different point on the Coulomb branch, involving two different masses. We 
extend this approach to single and double Regge limits of the five-gluon amplitude. In both cases, 
this alternative approach makes clear that log M5 should only have single logarithmic dependence 
on the kinematic variables because coUinear divergences are absent. A similar approach can be 
taken for some, but not necessarily all, of the n = 6 Regge limits considered in refs. [lOlllI]. 

The paper is organized as follows. In section 2 we study the four-gluon amplitude, presenting 
explicit results for all contributing four-loop integrals at the symmetric point s = t and arguing 
that, to all loops, the NLL contribution to the amplitude in the Regge limit is given by a small 
subset of all diagrams. In section 3 we turn our attention to the five-gluon amplitude at two 
loops, the evaluation of which leads to eq. (11. 5p . Section 4 contains a discussion of various Regge 
limits of Higgs-regulated amplitudes for n > 5, with emphasis on those features which differ from 
similar limits of dimensionally-regulated amplitudes. Section 5 summarizes our results, while 
various technical details can be found in three appendices. 



2 The four-point amplitude 

In ref. ^23j, it was suggested that the analog of the Bern-Dixon-Smirnov ansatz [12j for the planar 
four-point amplitude in Higgs regularization is 
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is the cusp anomalous dimension, and 

g,{a) = -C3a' + Oia'), 
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are analogs of functions appearing in the BDS ansatz in dimensional regularization [12], but 
need not be identical since they are scheme-dependent [23]. Overlapping soft and collinear IR 
divergences are responsible for the double logarithms in eq. (12. ip . The nontrivial content of 
eq. (12. ip is the statement about the finite terms; the IR singular terms of the amplitude are 
expected to obey eq. (12. ip on general field theory grounds (see refs. [i2lli3]). 

According to the assumption of dual conformal symmetry, the planar L-loop amplitude can 
be written as 

M(^) = J2 c(^) I (2.4) 



5 



where the sum runs over all extended dual conformal integrals /, with some coefficients c(/). 
In the case of the four-point amplitude, the coefficients are simply numbers. The set of loop 
integrals invariant under extended dual conformal symmetry is significantly smaller than that of 
generic loop integrals. 

At two and three loops, the assumption (12. 4p . together with the infrared consistency con- 
ditions, leads to a result in agreement with the exponential ansatz (12. ip . with the values (12. 2p . 
(1^ ancfl taOj 

g^^^ ^ 2.688870547851 ± 6.5 x 10"", cf ^ ^ -9.24826993 ± 9.6 x 10"^ . (2.5) 

If one assumes that these coefficients have transcendentalit}{§ 5 and 6 respectively, one finds for 
them the probable analytic valued 



Gi'' = - C2C3, 5f = -f Ca - 2Cs (2.6) 

using the PSLQ algorithm [i^H5]. 

The three-loop amplitude was computed [30j by assuming that 

Mf (s,t) = -^ \ha{s,t) + 2l3b{s,t)\ +is^t). (2.7) 

o L 

There are two additional dual conformal invariant three-loop integrals, /3c and I^d, with powers 
of in the numerator, but which have a finite m — > limit. These could in principle contribute 
to the three- loop amplitude if the coefficients multiplying them are nonzero. In ref. [30], we 
showed that even in this case the exponential ansatz remains valid, provided that the coefficients 
^Q^'* and cf^ are shifted accordingly. 

2.1 Four-loop four-point amplitude 

For the four-loop four-gluon amplitude, we will use the ansatz [32] 

^ {S, t) = ^ [ha{s, t) + 2/4fe(s, t) + 2/4,(5, t) + /4d(s, t) + 4/4e(s, t) + 2/4/(s, t) 

- 2hd, (s, t) - hf, (s, t)] +{s^t) (2.8) 

where the individual integrals, shown in fig.iH, are defined in dimensional regularization in ref. [32j . 
These integrals are all dual conformal invariant, and are straightforwardly rewritten in Higgs 



^In the course of the four-loop computation of the present paper we have improved the numerical accuracy of 
^Q^^ and quoted in ref. [30]. Here we display the results with improved numerical accuracy. 

®If we attribute a degree of transcendentality 0, 1, 1 and n, respectively to rational numbers, tt, log and Cn, and 
define the transcendentality of a product to be additive, then the i-loop amplitude is expected to have uniform 
transcendentality 2L. 

^We thank Lance Dixon for suggesting this value of ^q'^' to us. 
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Figure 1: The eight diagrams contributing to the four- loop four-point amphtude. We use the standard 
dual variable notation, labeling the external faces by xi through X4 and the internal faces by Xa through 
Xd- The former are related to the external momenta via pt = Xi — Xj+i (where i is understood mod 
4) while the latter are each integrated with the measure d^x/{m'^). Under each diagram is shown 
the numerator factor for the corresponding integral. To avoid clutter, we omit an overall factor of 
^13^24 ~ from each diagram (where Xab = ^a — Xb), and we do not label internal faces not appearing 
in numerator factors. As an illustrative example we demonstrate how to assemble all ingredients of the 
integral I^h in eq. (|2.9p . 



regularization, following refs. 
we have 



hb{s,t) 



. For example, using dual coordinates [Ij for convenience, 



(^2)4 



(xL + nl?){xl^^ + m^){xl^ + m?){xl^ + m^){xl^ + m^) 
1 



(2.9) 



Note that in 1^^ there is no +m? term in the loop-dependent numerator, since it connects two 
internal integration points. 

We have written down Mellin-Barnes representations for all eight integrals. This is easily 
done introducing the MB representation loop by loop [30lll6]. Interestingly, the dimensionality 
of most MB representations is one lower than the corresponding representation in dimensional 
regularization. At two and three loops, the opposite was the case. 
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To make contact with the exponential ansatz (12. ip . we will compute M^ls, t) at the symmetric 
point s = t. Defining x = w? /s = m? /t, the expression above becomes 



1 

16 L 



2ha{x)+Ah,{x)+Ah,{x) + 2hd{x) + %he{x) + Ahf{x)-Aha,{^)-If,{x) . (2.10) 



To evaluate these integrals, we proceed as in refs. [231 130]. The starting point is a multi- 
dimensional MB representation depending on the parameter x. For simplicity, consider a one- 
dimensional MB integral 

/3+ioo 

dzx' f{z) (2.11) 

— ioo 

where typically f{z) is a product of V functions, and /3 < 0. In principle, one could close the 
integration contour, say, on the right and obtain the answer as an infinite series arising from poles 
of the r functions in f{z). However, since we are only interested in the logo; terms as a; — 0, it 
is sufficient to deform the integration contour to positive values of Re(2;). The logarithms arise 
from taking residues at 2; = 0. In the case of multi-fold MB integrals, the above strategy can be 
iterated. We obtain expressions of the form YH^q h log^"^ x + C(log x) using the Mathematica 
code MBasymptotics [17]. In general, the coefficients bi still involve a (significantly lower) number 
of MB integrals, which we evaluate numerically using the code MB [JH]. Denoting L = logx, we 
find 

ha{x) = —L'' + —7t'L<'-—C3L'-—tt^L^ (2.12) 
^ 56 135 15^ 27 ^ ^ 

( 32 8 \ 

-yC2C3 - -C5 j + (-162.26621838645 ± 1.9 x 10"^°)^^ + 0{L) . 

/^.(a;) = iii^s - Ivr^L^ - -CaL^ + -vrU^ (2.13) 
^ ^ 5040 60 5^ 72 ^ ' 

+ (^C2C3 + ^C5^ + (525.46852427784 ± 9.9 x 10-''>)L^ + 0{L) . 

I^^(x) = ^IL^s - -CsL^ - ilvrU^ (2.14) 
^ ' 10080 3 270 ^ ' 

( 127 \ 

+ f -2C2C3 - — C5 j + (-128.86933736 ± 4.2 x 10"^)^^ + 0{L) . 

= ^L^ + —Ti^L^--C,L'>-—'K^L^ (2.15) 
^ ' 560 270 15 135 ^ ' 

(128 \ 
— ^(2(3 - 4OC5 J + (-710.51212126801 ± 1.4 x 10"^°)^^ + 0{L) . 

/^^(a;) = JZL^s + Ivr^L^ - 2^3^^ _ i^^4^4 ^2.16) 
^ ' 10080 60 3^ 2160 ^ ' 

23 33^ \ 

-7rC2C3 - -77(5 L"" + (222.7007725 ± 1.8 x IQ''^)L^ + 0{L) . 
6 12 / 
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'a;) = i^L^ - — TT^L^ - -C3^5 ^ -vr^L^ (2.17) 
^ ' 2520 135 5 ^ 30 ^ ' 

68C2C3 + ^Cs ) + (-117.32774717 ± 1.8 x 10-^)L2 + 0{L) . 



'4/ 



ImA^) = -^CsL' + ^n'L' (2.18) 
+ (^16(2(3 - yC5^ + (180.37203096920 ± 6.6 x 10"^°)^^ + 0{L) . 

/ (a;) = _ ii^2^6 _ 88^^5 ^ 11^4^4 ^2.19) 

■^'^ ^ 1260 135 15^ 15 ^ ' 

+ ^168(2(3 + -y^C5^ + (324.1906414642 ± 2.6 x 10"^)^^ + 0{L) . 

Summing up the contributions of the four-loop integrals computed above using eq. (12.101) . we 
obtain 

Mt\x) = ;^L^-C3i^^ + -^7rU^ + (6C2C3-9C5)i^' + (6.71603090 ±9.1 x10-')L2 + O(L). (2.20) 
24 60 

The coefficients of L"^ and in the expressions above were obtained numerically, so we cannot 
distinguish between rational and transcendental numbers. Nevertheless, motivated by the ex- 
pectation that the result should have uniform transcendentality, we have replaced the numerical 
values by their probable analytical equivalents. Specifically, we used the Mathematica implemen- 
tation of the PSLQ algorithm to identify linear combinations of numbers with the correct degree 
of transcendentality that agree with our results within the numerical accuracy. It goes without 
saying that this does not constitute a proof that these expressions are necessarily correct. A 
guess for the numerical coefficient of the term is — ICe + SCf- 
Exponentiation of the IR logarithms requires that 

1 ... 1 



Mf(x) = - C3l^^ + ^vrU^ + (4C2C3 - 2^J^) 



- - + 2Ci) + Oil) (2.21) 

where we have used the values of 7(a), Go{a), and £4(0) given in eqs. (12.21) and (12.31) . There is 
complete agreement between eqs. fl2.20p and fl2.2ip . Comparing the coefficients of L^, we confirm 
the value ^q^^ = ICs — C2C3 obtained at three loops by assuming that I^^ does not contribute. 
Comparing the coefficients of L^, we confirm the value cf^ = — ^Ce ~"2^| obtained at three loops 
by assuming that I^c does not contribute, provided that the four-loop cusp anomalous dimension 
is given by 

7(4) = -117.1788222 ± 3.7 x 10^^ ^ -4(2 - 24(2(4 - 4C| - 50(6 (2.22) 

in perfect agreement with the result (also numerical) found in refs. [321E3], and in agreement 
with the spin chain prediction from ref. [49]. It is noteworthy that our result (12.221) improves 
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the numerical precision by two orders of magnitude compared to ref. [33] and by five orders of 
magnitude compared to ref. [32j^. 

Similar to the situation at three loops, there are 10 four-loop integrals (in addition to those 
shown in eq. f l2.8p ) that could in principle contribute to the four-loop amplitude, based solely on 
requiring dual conformal invariance |50]. As we have seen, however, our results are completely 
consistent with the absence of these additional integrals at three and four loops. 

It is amusing that the coefficients appearing in the results for the individual integrals are 
rather complicated (e.g. the coefficients of the terms), yet they sum up to give the very 
simple result f l2.20p . as required by infrared consistency. This suggests there may exist a better 
organization of the calculation that avoids the complexity of the intermediate results. 



2.2 Regge limit of the four-point function 

The four-point amplitude (12. ip can be rewritten as 

1 ~ vr^ ^ 

logM4 = --7(a) (log t;) (log m) + g^ia) {\ogu + logf ) + y7(a) + C4(a) + C(m^) 



(2.23) 



where u = m? / s and v = m^/t. Note that (logw)^ and (logf)^ terms are absent, which is related 
to the Regge-exactness of the four-point amplitude. This can alternatively be interpreted as 
follows [30j . The Higgs-regulated amplitude can be computed at a different point in moduh space, 
where the diagrams contain internal lines with different masses m and M along the periphery, 
and external lines with mass M — m. By dual conformal invariance, the amplitude only depends 
on M = m^/s and v = M'^/t. The absence of (logm^)^ terms in the amplitude can be understood 
as the absence of collinear divergences when massive (M) particles scatter by exchanging lighter 
(m) particles (with the mass of the lighter particles serving as an IR regulator). 

We may exponentiate eq. (12.231) . and expand the result in powers of the coupling a and of 
the Regge logarithms logw. At L-loop order, the leading logarithm is log^w, NLL is log^ 
etc. To LL and NLL order, the L-loop amplitude in the Regge limit is given by 
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(2.24) 



where we have used eqs. ( 12. 2p and ( 12. 3p . In contrast to dimensional regularization, in Higgs 
regularization there is a single diagram, the vertical ladder Iiaiv, u) (see fig.E]), that contribute^ 
to the LL term of the L-loop amplitude (I2.24p . In the LL limit, the vertical ladder factorizes 
into a product of (two-dimensional) bubble integrals (again see fig. [2]). Moreover, the LL, NLL, 
and NNLL contributions of the vertical ladder diagram were computed (cf. eq. (4.16) in ref. |30] ) 



^°The value (|2.22l) corresponds to r = — 1.99999892±6.3 x 10 ^ in the parameterization used in those references, 
"'^^in the Regge (b) hmit (see footnote [6]) . 
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Figure 2: Factorization of the leading- log and next-to-leading-log contributions to the Regge limit 
s ^ t of the L-loop vertical ladder integral lLa{v,u) into simpler integrals. Factorization of the NLL 
contribution of the vertical ladder integral with H-shaped insertion Ilh- The dotted line indicates a 
loop-momentum-dependent numerator . 

using the method of ref. [51]. Subtracting the vertical ladder contribution from the prediction 
( I2.24P for the full amplitude, we obtain 

(2.25) 

From this, one sees that contributions from diagrams other than the vertical ladder are required 
at NLL order to obtain the expected amplitude in the Regge limit. 

Rules for evaluating the leading logs behavior of multiloop integrals were summarized in 
refs. [511152]. One begins by identifying paths through the graph, which, when contracted to a 
point, split the diagram into two parts with a single vertex in common, and with pi and on one 
side and p2 and ps on the other (for example, each of the rungs of the vertical ladder diagram). 
Paths of minimal length are called "d-lines" [53] or "t-paths" [51] • A scalar diagram containing 
m d-lines of length n goes as log™~^ s/s" as s — ?■ oo. For example, the L-loop vertical ladder 
diagram contains (L + 1) d-lines of length one, so the vertical ladder integral (multiplied by st^ 
to make it dual conformal invariant) goes as log^s, or equivalently log'^w, as discussed above. 
All other L-loop diagrams contain at most (L — 1) d-lines and hence prima facie give at most 
a log'^~^M contribution. For example, the four-loop diagrams hciv,u) and hdiv^u) (see fig. [1]) 
each contain two d-lines of length one and prima facie go as logw. (These diagrams also contain 
two paths of length two, which are not minimal and therefore do not contribute.) 

The MB calculation summarized in appendix |A1 however, shows that both of these diagrams 
go as log^ u, two powers higher than expected. This is because the d-line rules given above 
only apply to scalar diagrams with no non-trivial (i.e., loop- momentum-dependent) numerator 
factors. The presence of numerator factors, which are required for the loop integrations next to 
the H-shaped insertion to have the correct dual conformal weight, can increase the leading power 
of log u of the diagram. The reader may ascertain from the results of appendix |A] that diagrams 
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with no nontrivial numerator factors have the logM dependence predicted by the d-hne rules, 
whereas those with numerator factors can have a stronger logw dependence. 

We show in appendix |B] that when an H-shaped insertion in a vertical ladder diagram is 
accompanied by a numerator factor, the two lines of length two constituting the sides of the H 
are effectively promoted to length one, increasing the d-line count by two. In particular, an L- 
loop vertical ladder diagram with a single H-shaped insertion I^h (see fig. [2]), which prima facie 
would go as log^~^ M, actually goes as log'^"^ u due to its numerator factors, and thus contributes 
to the amplitude at NLL order. A calculation in appendix |B] further shows that, subject to 
reasonable assumptions, the leading log contribution of the integral Ilh factorizes as 

/LH = |^^log^-'wxK(t;)^-2^ir'(t;) + 0(log^-2w), (2.26) 

where K{v) and K'{v) correspond to the two-dimensional bubble and two-loop bubble diagrams 
shown in fig. |2] (see ref. [30] for further discussion). Taking v small, we have 

K{v) = -2logv + 0{v), 

K'{v) = -^log^v -^n^logv + 0{v) . (2.27) 

Note that eq. fl2.26p implies that the position where the H-shaped insertion is made into the 
vertical ladder integral is unimportant at NLL order, i.e. all such integrals give the same NLL 
contribution. 

At two loops, Ilh is just the horizontal ladder l2a{u,v). At three loops, Ilh is the tennis 
court diagram I^biu, v). In ref. jSOj, it was shown that the leading contribution of these diagrams 
is precisely given by eq. f l2.26p . At four loops, the three Ilh diagrams are given by /4c(f,M), 
its flipped version, and l4d{v,u). The results for these integrals found in appendix lAl are also 
in agreement with eq. f l2.26p and, moreover, all other four-loop diagrams contribute at most to 
NNLL. Finally, we have verified eq. (I2.26P at five loops as well using the integrals given in ref. [5]. 

At L-loop order there are (L — 1) vertical ladder diagrams containing a single if-shaped 
insertion, so that the total contribution to the amplitude of the vertical ladders with one H- 
shaped insertion is (multiplying the contributions by (—1/2)^) 

' {L - 1)Ilh = ^^^T^ log^"' « [log"-^' ^' + log"-"' ^ + 0{v)] + 0(log^-^ u) (2.28) 

which precisely matches the result fl2.25p expected from the exponential ansatz. 

In summary, we have shown that in Higgs regularization, the NLL contribution to the four- 
gluon amplitude in the Regge limit is given by a small set of diagrams: the vertical ladders and 
the vertical ladders with one H-shaped insertion. We confirmed this through five loops by direct 
evaluation of the integrals, and we gave an argument that this holds to all loop orders. 



3 The five-point amplitude 

In this section, we direct our attention to the two-loop n = 5 point amplitude of A/" = 4 SYM 
theory, which is of interest for several reasons. It serves as a further instructive application 
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P3 



Figure 3: Dual conformal scalar integrals contributing to the five-particle amplitude at one and 
two loops. The dashed line indicates that the integral contains the loop-momentum-dependent 
numerator factor {q + p^Y + m^. 

of the Higgs mechanism to regulate infrared divergences, and confirms the universality of the 
exponential structure of IR singularities of n-point amplitudes in massless gauge theories (which 
is usually studied in dimensional regularization; for example, see refs. [i2| H3]). We also establish 
an iterative relation at two loops, which is the exact analog of the five-point iterative relation 
in dimensional regularization ^55] EH]- This allows us to write an all-loop ansatz for n-point 
amplitudes in Higgs regularization, analogous to the BDS ansatz in dimensional regularization, 
whose content is that the IR-finite part of the amplitudes also exponentiates. Of course the 
separation between IR-divergent and IR-finite terms is not unique since one could always add a 
constant to one while subtracting the same constant from the other. Knowledge of the four-point 
amplitude alone does not give enough information to resolve this ambiguity in a natural way, but 
after computing the five-point amplitude, we will be able to determine a unique way of writing 
the universal IR-divergent part of the Higgs-regulated scattering amplitude for any n. 

As in ref. [30] we evade a Feynman diagram calculation by beginning with the ansatz that 
Higgs-regulated five-point loop amplitudes (normalized, as usual, by dividing by the correspond- 
ing tree amplitudes) can be expressed as linear combinations of all possible dual conformally 
invariant scalar integrals. At one loop there is a unique integral I^'^\si) (see fig. [3]) that can 
appear in the ansatz 

^s'^ = -\H ^i^54'^(3.) + Oim'), (3.1) 

cyclic 

where +(9(m^) stands for potential parity-odd terms (see below), 

Si = {Pi+Pi+if, i = l,...,5 (3.2) 

(with pe = pi), and the sum in eq. (13.11) runs over the five cyclic permutations of the external 
momenta Pi. The small expansion of I^^\si) is given in appendix [Cl After summing over 
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cyclic permutations, the one-loop amplitude simplifies to 



m: 



(1) 



E 

1=1 



1 1 2 f Si 

4 Vm"' 



(3.3) 



where the corresponding finite remainder F^^^ is given b}@ 



1 ^ 



1=1 



log 



log 



-51+2 



TC 



(3.4) 



with Si+5 = Sj. 

At two loops there are two dual conformal invariant scalar integrals that contribute to the 



amplitude: the double box I^^^"'{si) and the pentagon-box I^'''''{si) with a numerator factor 
involving the pentagon loop momentum (see fig. [3]). The coefficients of these integrals in the 



{2)b, 



amplitude Mg 

ansatz 



(2), 



M, 



(2) 



Si) are determined by the consistency of infrared singularities, leading to the 
1 



cyclic 



(2)a, 



I 2 r(2)a/ 



+ S2S3S5/, 



(2)6, 



+ 0{m' 



(3.5) 



where, as in the four-point calculation, the relative numerical coefficients of each diagram are 

precisely the same as in dimensional regularization. 

Out of abundance of caution we have included +0(m^) in both eqs. (13. 5p and (13. ip in order to 

(2) 

encapsulate possible parity-odd terms. In dimensional regularization, M5 has a non-vanishing 
parity-odd contribution starting at 0{e~^). The parity-odd part of logMs, however, is 0{e) at 
one- and two-loop order [551457] due to the non-trivial cancellation of parity-odd terms in M5 

and — ^(Mg^'')^, which also has a contribution starting at 0{e~^). The vanishing of the parity- 

(2) 

odd contribution to the Higgs-regulated amplitude Mg , together with the fact that terms that 
vanish as m — )• are not required to compute logMg in Higgs regularization, implies that the 
parity-odd contribution to logMg is 0{m'^) at one and two loops. 



The small m expansions of Jg 
over cyclic permutations, we find 



(2)a, 



Si) and /5^'"'(sj) are given in appendix O After summing 



i=l 



Using the identity 

5 



Ei°f 



i=l 



C2 , 2 f Si\ C3. f Si\ (2 2 



4 



Si 



io<^)-|:[2io.e)iog(^)-2io, 



C2 



Si+l/ 4 \Si+2 



2 / log 



(3.6) 



(3.7) 



Recently it has been shown [581159) that, when expressed in momentum-twistor variables, the expression in 
eqs. (I3.3P and (I3.4p computes the volume of a 4-simplex in AdS^ with 5 (regulated) points on the boundary. 
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we may rewrite this as 



M, 



(2) 



- M 



(1) 



i=l ^ 



■^log^f^UllogfiL 



where c^^ = |C4 and F^^'{si) is given in eq. fl3.4p . Equation (I3.8p is analogous to the two-loop 



five-point iterative relation of refs. 

It is instructive to rewrite the one- and two-loop four-point amplitudes in a form similar to 
eqs. fl3.3p and fl3.8p . Using eq. ( I2.ip . with s = si = and t = S2 = S4 for four-particle kinematics, 
we have 



(1) 



M, 



(2) 



I {Mr 



i=l 
4 

E 

i=l 



-7 log 



4 



C2 , 2 
J log 



+ ttIos 



- (,2 -^4 + C4 



where c 



;:(2) 



[(^4 and 



i=l 



log^ 



By re-expressing the constants in eqs. (13. 8p and fl3.10p as 



;;(2) 



{2n - 5) 



C4 



we may combine eqs. (13. 8 p and (I3.10p into 



n 



i=l 



C2 



^log^P^ +^logM7 +^ 



-C2F«-^C4 + OK) 



Note that the common infrared divergence term 

n 



4 



resembles the form expected from dimensional regularization |l2lll3l[60] : 

- [C2 + C36 + C4e']M«(2e). 



(3.9) 
(3.10) 



(3.11) 



(3.12) 



n = 4, 5. 
(3.13) 

(3.14) 
(3.15) 



It would be interesting to understand more precisely the relation between these two forms of the 
infrared divergences, perhaps along the lines of ref. |61jl^ 
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We are grateful to S. Moch for discussion and correspondence on this question. 
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Armed with the above ingredients we are now in a position to pose the n-point generahzation 
of eq. (El]) as 



logM„ = ^ 
where 



i=l 



■4^iog=(4 

lb Vm"^ 



^ log )+/(«) 



1^ 
4 



+ -7(a) F^i) + 7^„ + C'(a) + ^(m^ 



/(«) 



C{a) 



(3.16) 
(3.17) 



and the remainder function TZn vanishes for n = 4 and n = 5. 

We can then proceed by defining the Higgs regularization analogs of the IR-finite functions 
Fn introduced in ref. [T2] by subtracting the universal infrared singularities from the logarithm 
of the amplitude 



Fn= lim logM„-V 

m2-i.O > 



-^iog=(4 

16 Vm"^ 



— — log(— + / a 

2 \m'^ 



(3.18) 



The precise forms of the one-loop functions (13.41) and (13.111) and the two-loop functions 

35C4" 



.(2) 



.(2) 



4 

E 

i=l 


C2 , 2 

log 

8 ^ 


5 

E 

i=l 


C2 , 2 
-^log 



Si+l 



16 

+ ^ I0-2 



4 



Si+2 



3C4 

2 



(3.19) 



differ (by additive constants) from the corresponding expressions in dimensional regularization, 
which may be found in ref. [12]. Through two loops, however, the n = 4 and n = 5 IR-finite 
functions satisfy a similar iterative relation 



F^^) + C{a) + 0{a!' 



= 4, 5 . 



(3.20) 



While the constants Qo{a) and c„(a) defined above take different values in Higgs and in dimen- 
sional regularization, it is fascinating to note that the two-loop value of C{a) in eq. (I3.16P is 
identical to the corresponding value in dimensional regularization |i60j. Perhaps this is a coin- 
cidence, or perhaps it has a deeper explanation, especially in light of the fact that the same 
value appears yet again in the finite part of the two-loop lightlike polygon Wilson loop after 
appropriate subtraction of U V divergences [52] ■ 

For n > 6, we expect that the "remainder" function TZn in eq. ( I3.16P is non-trivial, just 
as in the corresponding formula in dimensional regularization [T5l[T6]. However it is natural to 
expect TZn to take the same value in Higgs regularization as it does in dimensional regularization. 
This is because the remainder function is an infrared-finite, dual conformally invariant quantity 
(as required by the dual conformal Ward identity), so it constitutes a good "observable" of 
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SYM theory. The 6-particle remainder function at two loops was first computed numerically 
in ref. [TJ], where agreement with the corresponding remainder function for lightlike hexagon 
Wilson loops [16l[63] was established, and an analytic expression has been given more recently 
in refs. [IHM]. 



4 Regge limits 

In this section we discuss a number of features of the Regge limits of Higgs-regulated amplitudes 
for n > 5. As emphasized in footnote [6] such limits may be taken in two different orders: limits (a) 
where all mf are first taken to be much smaller than all kinematical invariants, and subsequently 
a Regge limit is taken, and limits (b) where the Regge limits are taken first with various fixed 
kinematic invariants and fixed masses , which are subsequently taken to be much smaller than 
the fixed kinematic invariants. The Regge behavior of individual diagrams can depend on the 
order, (a) or (b), in which these limits are taken. 



4.1 The five-point amplitude 

For the purposes of discussing the Regge limits of the five-gluon amplitude, we adopt the following 
parameterization of the kinematical invariants |1D] 

s = (p^+p2)^ ^l = (P2+P3)^ ■sl = (p3+P4)^ 

S2 = (Pi+PbY, ^2 = (P5+Pl)^ S = SiS2/k. (4.1) 

We will consider two different limits of these invariants, as defined in ref. |40] : 

Single Regge limit : s ^ oo, Si — i- oo, k, S2, ti, t2 fixed , (4.2) 

Double Regge limit : s — ?■ 00, si — )■ 00, S2 — > 00 K,ti,t2 fixed. (4.3) 

In the following two sections, we will consider the Regge limits (a) and (b) for the five-point 
amplitude. We will see that, at least through two loops, the amplitude is the same in both 
limits. 



4.2 Regge (a) limits for n = 5 

If we wish to take Regge limit (a) of the five-gluon amplitude, we can start with the conjectured 
ansatz f l3.16p . in which the small limit has already been taken. In terms of the parameters 
( 14. ip . eq. f l3.16p for n = 5 takes the form 



log Ms = w(ti)log(^ 



+ 



16 



7(a) 



log 



log^ 



i0(t2 
tl 

\- log 



log 
-log^ 



+ 2 log (^) log 



tlt2 



- log 



K 



(4.4) 
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where 

1 , . . / t 



uj{t) = a{t) - 1 = --7(a) log(— ) - ^o(a) (4.5) 
is the same trajectory as in the four-point function (11.31) . Equation (14. 4p is equivalent to [T^HO] 



M.H^)""''(4)*'f('..«^.M. (4.6) 



which exhibits the expected factorization. From eq. (14. 6 p it is straightforward to take the Regge 
hmit to obtain the following expressions, separating the Regge behavior from the fixed term: 



Single Regge (a) limit : M5 — > — - — - F(ti,t2, k, 



S2 \'^(*2) 



(4.7) 



^ M7 \F(ti,t2,K)] . (4.8) 

Equations (14. 7p and (14. 8 p have the same form as in ref . |10] , but with a different value for the 
constant term in the trajectory function (14.50 . and a different F{ti,t2,K) due to the difference 
between eq. (14. 4p and the analogous BDS result. 

4.3 Regge (b) limits for n = 5 

To obtain the Regge (b) limits of the five-point amplitude, one must start with scalar integrals 
with finite m, take the Regge limit first, and only afterwards take m small. We have evaluated 
the five-point one- and two-loop diagrams in this way and have obtained results identical to 
eq. (14. 7p in the single Regge limit and to eq. fl4.8D in the double Regge limit. Hence, at least to 
two-loop order, the amplitude is independent! ^^1 of the order in which the Regge limit is taken 
(even though the individual diagrams are not). 

In ref. [30], an alternative approach to the Regge (b) limits was applied to explain the absence 
of double logarithms in the Regge limit of the four-gluon amplitude. We can do the same for 
various Regge (b) limits of five-gluon amplitudes. Consider the theory at a different point on the 
Coulomb branch, where the scalar diagrams have internal lines with variable masses along the 
periphery (and vanishing masses in the interior). Let rrii be the mass of the line(s) connecting 
and Pi. The masses of the external lines are given by = — (m^ — rrij+i)^. Due to dual 
conformal symmetry, the amplitude depends only on the dual conformal invariants 

mjmj_|_2 

Ui,i+2 = 7 ■ r^—T r^. (4.9) 

[Pi + Pi+i)^ + [mi - mi+2Y 

By making various choices for mj, we can reproduce the single and double Regge limits. 



In an "unphysical" Regge limit s ^ 00 with si, 52,^1,^2 fixed, described on p. 177 of ref. [5T], the one-loop 
amplitude does depend on the order of limits, yielding log^ s dependence in the Regge (a) limit, as can be seen 
from eq. (|4.4p . but logs dependence in the Regge (b) limit. 
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Figure 4: Two-mass configurations illustrating the (i) single and (ii) double Regge (b) limits for the five- 
point amplitude. Fat and thin lines along the periphery denote particles of mass M and m respectively. 
Fat and thin exterior lines denote particles of mass M — m and respectively. 



Single Regge (b) limit: 

By making the choice = m, and rrii = M for i ^ 3, we have 

Mm _ _ Mm _ _ 

s + {M-mr s, + {M-mr ~ W ~ ~ ' 

(4.10) 

Then taking the hmit m <^ M yields misjMss <C M24, 'U41, ^52, which is equivalent to the single 
Regge limit (14. 2p . The resulting diagrams (e.g., see fig. Ilt^i)) cannot have coUinear divergences 
because the massless external lines never connect to the light mass internal lines. The amplitude 
therefore has at most a simple logm IR divergence, which corresponds to logwas (or logwis) and 
therefore to a simple logsi dependence in the single Regge limit, in agreement with eq. (14. 4p . 

Double Regge (b) limit: 

By choosing mi = m^ = m, and m2 = m^ = m^ = M, we have 

Mm Mm M^ 

= - ^^^ = 77' Si + {M-mr = s, + {M-mr = ~ ' 

(4.11) 

Taking the limit m <^ M yields -Uia ^ M35,M4i ^ M24,W52, which is equivalent to the double 
Regge limit (14. 3p . Again, none of the diagrams that contribute to the amplitude in this limit 
has collinear divergences (e.g., see fig. IH^ii)), hence the amplitude has at most a simple logm 
IR divergence, which corresponds to logwas or logM4i, and therefore to a simple logsi or logS2 
dependence in the double Regge limit, in agreement with eq. (14. 4p . 



4.4 Regge (a) limits for n > 6 

To consider various Regge (a) limits of n > 6 amplitudes, we may start with the conjectured 
ansatz (I3.16p . in which the small m^ limit has already been taken. 

For n = A and n = 5 we have seen that Fn \ the IR-finite part of the one-loop amplitude, has 
the same dependence on kinematic variables (up to an additive constant) as the corresponding 
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function in dimensional regularization. Any difference in the form of F„ for n > Q between 
Higgs-regulated amplitudes and dimensionally-regulated amplitudes should be dual conformal 
invariant, and therefore a function of dual cross-ratios. Any such function remains finite in 
Euclidean Regge limits |1D]. Similarly, the remainder function in either Higgs regularization 
or in dimensional regularization remains finite in Euclidean Regge limits, as discussed in sec. 7 
of ref. [lO]. Therefore, the Regge (a) limits of Higgs-regulated amplitudes are equivalent to the 
Euclidean Regge limits of dimensionally-regulated amplitudes discussed in refs. [IQIIIT] , but with 
the Regge trajectories given by eq. (14. 5p . Similarly the Regge vertex functions will be analogous 
to those in ref. |1D], but will differ in detail due to the difference in the IR- regulator scheme. 

It is possible that Tin could contribute to Regge limits in the physical region, as these may 
involve contributions which do not remain finite. However, it is known that there are difficulties 
in continuing M„ [n > 6) from the Euclidean to the physical region [TllHTl|6l]. These same 
difficulties would be present for Higgs-regulated amplitudes, so we do not consider these latter 
limits further. 

In conclusion, there are no important differences between the Regge (a) limits of Higgs- 
regulated vs. dimensionally-regulated n-gluon amplitudes. 

4.5 Regge limit (b) for n = 6 

We have not attempted to compute the amplitudes for > 6 in any of the Regge (b) limits 
defined in pojl^ to ascertain whether they are equivalent to the Regge (a) limits discussed in 
the previous subsection. 

However, we would like to point out that the analysis of some of these limits could be facili- 
tated by considering Higgs-regulated amplitudes at various points on the Coulomb branch. Let 
rrii be the mass of the lines on the periphery of a diagram connecting external lines and pt, 
with vanishing masses in the interior. We consider points on the Coulomb branch involving just 
two distinct mass assignments m and M (with m <^ M). 

There are four inequivalent configurations (for six-point amplitudes) that yield at most single 
logarithmic dependence on the IR-regulator mass m. (These are configurations where no two 
adjacent rrii are equal to m.) We can classify these configurations according to the number of 
lines with the small mass m. 

One line with small mass. 

If we set mi = m and 1712 through mg equal to M, then none of the massless external lines are 
attached to the line(s) on the periphery with mass m. Hence there are no collinear divergences, 
and the amplitude only goes as logm. 

Two lines with small mass. 

There are two possibilities involving two lines with small mass. The first one consists in 
setting nil = "^3 = with the remaining masses given by M. The second consist in setting 
mi = 7714 = m, with the remaining masses given by M. The latter corresponds to the Mueller- 
Regge limit discussed in section 7.4 of ref. [H]. 
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Three lines with small mass. 

If the external masses alternate between m and M (i.e. mi = = = m and m2 = = 
rriQ = M), one obtains the poly-Regge limit of section 7.3 of ref. |41j . 

No two- mass set-up has been found to give the single Regge limit or helicity pole limit. 

5 Summary 

In this paper we have continued the program of computing higher-loop A/" = 4 SYM planar n- 
gluon amplitudes and testing various conjectures using the Higgs regulator scheme proposed in 
ref. [23] and further developed in ref. [30] • Specifically, we have extended the analysis to four loops 
for the four-gluon amplitude, and to two loops for the five-gluon amplitude, using Mellin-Barnes 
techniques to evaluate the integrals. We have assumed that only scalar diagrams invariant under 
extended dual conformal symmetry contribute to the amplitudes, and with the same numerical 
coefficients as in dimensional regularization. Although one has no a priori guarantee that the 
set of diagrams contributing to the amplitude in one IR-regulator scheme coincides with the set 
in another scheme, all the results we obtained using this assumption are consistent with the 
universal IR-divergence structure of massless gauge theories and also with the conjectured all- 
loop ansatz (11.71) for the IR-finite part. For all cases considered, we have verified that the IR- finite 
parts of the logarithm of the amplitudes have the same dependence on kinematic variables as the 
corresponding functions in dimensionally-regulated amplitudes (up to overall additive constants, 
which we determine). 

We have also extended the study of the Regge behavior of Higgs-regulated amplitudes which 
was begun in ref. [SD]. The Regge (a) limits (in which the masses are first taken much smaller 
than all kinematic invariants, and then the Regge limit of kinematic variables is applied) can be 
understood by simply taking the kinematic limits of the ansatz (11. 7p . in which 0{m?) terms are 
already neglected. Various (a) type limits are discussed for n > 5 which give essentially the same 
results as those from the dimensionally-regulated BDS ansatz, but with different expressions for 
the gluon trajectory and Regge vertices resulting from the different regulator scheme. 

To study Regge (b) limits (in which the kinematical limits are first taken with fixed regulator 
masses, which are subsequently taken to be much smaller than the fixed kinematic invariants) 
one must evaluate the Regge limits of the individual diagrams contributing to the amplitudes. In 
the Regge (b) limit, certain classes of diagrams are dominant, whereas in the Regge (a) limit no 
single class of diagrams dominates. In ref. [30j, it was shown that the leading-log approximation 
of the four-gluon amplitude is dominated to all loop orders by the sum of vertical ladder diagrams 
only. In this paper, we showed that the next-to-leading-log approximation depends only on the 
vertical ladder diagrams, together with the class of vertical ladder diagrams with a single H- 
shaped insertion. 

One way of analyzing Regge (b) limits of diagrams is to go to a different point on the Coulomb 
branch involving several different masses. Examples are given for several Regge limits of n = 5 
and n = Q amplitudes using two different masses, although we have not been able to obtain all 
Regge (b) limits using two-mass configurations on the Coulomb branch. Although the Regge (a) 
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and Regge (b) limits of individual diagrams differ, we have found that the full amplitudes are 
independent of the order of limits in the cases that we have considered. 

The results of this paper show that the Higgs regulator for planar A/" = 4 SYM amplitudes 
continues to exhibit a number of practical and conceptual advantages compared to other reg- 
ulators, the first signs of which were observed in ref. [23l[30]. On the practical side, the Higgs 
regulated multi-loop integrals we have encountered so far have proven quite a bit simpler to 
evaluate than their counterparts in dimensional regularization. One consequence of this is that 
we have been able to compute the four-loop cusp anomalous dimension with numerical precision 
five orders of magnitude greater than ref. [32| and two orders of magnitude greater than ref. [33] . 
The crucial conceptual advantage of the Higgs regulator is that it preserves the remarkable (ex- 
tended) dual conformal symmetry which has recently played such an important role in unlocking 
the hidden structure of SYM amplitudes. While dimensional regularization seems completely at 
odds with modern (inherently four- dimensional) twistor-space methods, it is greatly encouraging 
that the Higgs regulator can be very naturally implemented in momentum twistor space, as seen 
for example in the beautiful recent results of ref. [58|l59]. which we suspect are just the tip of the 
iceberg. 
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A Regge limits of four-loop four-point integrals 

In this appendix, wc list the results for the Regge (b) limit of the each of the four-loop in- 
tegrals that contribute to the four-loop amplitude. For compactness, we employ the notation 
{ai, . . . , a„} = El=i «m log"-"^ v + 0{v). 

hais, t) = log ^ 0, ^tt', 0, ^7^^ 0, 252.30 . . . , const} + 0(log° u) . (A. 1) 

74a(i,s)=log^«{|o,0,0} (A.2) 
+ Iog3«{-|o,-^7r2, -8(3,0,0} 

+ ^°^^-^f'°'^-^'y^-f.M83.07..., 11.55...} 

X 352 736 68 608 
+ log« {-^, 0, - 4^7r2, --C3, -580.84 . . . , -1536.93 . . . , const} + 0(log° u). 

h,{s, t) = V u 0, ^7^^ 0, |7^^ 0, 0} (A.3) 

+ "^"1^' °' "i^""'' "Y "li""'' "^^-^^ ■ ■ ■ ' ■ ■ ■ ' + ■ 

/4fe(t, s) = logM {|-, 0, ^TT^, 0, |-7^^ 0, 65.11 . . . , const} + 0(log° u) . (A.4) 
63 45 z7 

16 8 32 

/4c(s, = log {— , 0, -TT^, 0, — 7^^ 0, 105.80 . . . , const} + C»(log° u) . (A.5) 
oio 45 io5 

/4e(t,s) =log'-^ u{^,0,^7r2, 0,0,0} (A.6) 

, 2 r 121 56 p 32, 32 , 
+ V u {-—, 0, — ^7^^ -y C3, - Y^7^^ -56.52 . . . , 0} 

,1963 539 5 62, 1478 . . n ^ 

+ logtz {-^, 0, — vr^, -Cs, ^7^^ 622.32 . . . , 1619.96 . . . , const} + 0(log° u) . 

hd(s, i) = log « {^, 0, ^vr^, 0, ^7^^ 0, 252.30 . . . , const} + 0(log° u) . (A.7) 

hdit, s) = log' u{^,0, ^TT^, 0, 0, 0} (A.8) 

+ V u {-|, 0, -|7^^ -SCa, -f ^^ -49.77 . . . , 0} 

+ \ogu{^, 0, 127r2, ^Cs, -191.36 . . . , 863.63 . . . , const} + 0(log° u) . 

8 16 56 

J4e(s, t) = log w { 0, — TT^, 0, — 7^^ 0, 130.22 . . . , const} + C»(log° u) . (A.9) 

he{t, s) = log' ^ 0' 0' 0' 0} (A-10) 

223 143 16 22 
+ log« {-^, 0, -^tt', --Cs, -y 7^^ -49.77 . . . , -258.40 . . . , const} + 0(log° u) . 
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hf{s, t) = log« 0, ^TT^, -^C3, -39.93 . . . , 41.35 . . . , const} + C»(log° u) . (A.ll) 

hf{t, s) = V w 0, yTT^, 0, ^7^^ 0, 0} (A.12) 

484 2fiS 734 
+ log^ {- 3J^' 0' -12C3, - Y3^7r^ -103.20 . . . , -907.47 . . . , const} + O(log° u) . 

4 2 

^4d2(s, t) = logM {0, 0, 0, --Cs, T77^^ -39.93 . . . , 41.34 . . . , const} + C(log° u) . (A.13) 

74d,(i, s) = log {0, 0, 0, -^Cs, 4^"' -9-84 ■ ■ ■ , 28.48 . . . , const} + e»(log° u) . (A.14) 

3 15 

hf,{s, t) = log' u{^,0, ^7r\ 0, ^7^^ 0, 0} (A.15) 

124 236 136 
+ log^ {-^' 0' "15"'''' ^ ZT""'' ■ ■ ■ ' -935.95 . . . , const} + 0(log° u) . 

hf,{t,s)=hf,{s,t). (A. 16) 
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B Regge limit of the L-loop ladder with //-insertion 



In this appendix, we evaluate the leading log contribution (in the limit of large s) of the dual 
conformal invariant L-loop vertical ladder diagram with one H-shaped insertion, Ilh, described 
in sec. 12.21 

We begin by considering an {L — l)-loop vertical ladder diagram, with the external regions 
labeled by Xi through X4 (where s = x^g and t = x^^), and the loops labeled by Xj, with 
i = 5, - ■ ■ , L + 4 (see fig. Et^i)). We replace the jth loop with a horizontal double loop, converting 
it into an L-loop diagram, as shown in fig. |5]^ii). By dual conformal invariance, this diagram 
must be accompanied by a factor of 

We now perform the integration over the double box. Its Feynman parameterization is 



'double box 



f ^ ^ ^ ^ ^ 2A(5(ao + ai + "2 + 51 + ^2-1) „^ 
/ daodaida2dpidp2djidj2 tt^ ^ — tt^ (B.2j 



where 



D = D24t + Dj^ij^iPj^ij^i + -Dj-1,2-Pj-1,2 + -Dj+i_2-Pj+l,2 + Dj-l,4Pj-lA + Dj+l,4Pj+l,4 

A = aoai + aoa2 + aia2 + ao5i + ai{5i + 62) + 0262 + 5i52 
a = ao + ^2 

S^ = Pi + ^i (B.3) 

with Pij = xjj = {xi — XjY . 

The coefficients of Pit in eq. ( 1B.3|) are given by 

D24 = aoaia2, ^j-ij+i = ("2 + 52)/3i7i + ("o + 5i)/3272 + ai(/3i + /32)(7i + 72), 
-Dj-1,2 = ao(ai/3i + "2/31 + ai/32 + 52/3i), = ao(ai7i + "271 + «i72 + (^271), 

L'j_i,4 = a2(ai/32 + ao/32 + + 5i/32), ^i+i,4 = «2(ai72 + «o72 + «i7i + '^172) • 

(B.4) 

If the double box is inserted at one of the ends of the vertical ladder, the expression for a in 
eq. (]B.3p will contain additional terms f3i + (32 or 71+72, but the results below will be unaffected 
by this change. 

We assume, following chapter 8 of ref. [65], that the asymptotic behavior of multiloop integrals 
is dominated by a region in the space of the loop momenta that can be described by a set of 
nested inequalities. In our specific case, those inequalities imply that the leading s — >■ 00 behavior 
of Ilh comes from the region where Pj-ij^i is large (i.e. a nonzero fraction of s) with the other 
Pik appearing in eq. (IB.Sp remaining finite. This in turn implies that the main contribution to 



the subintegral (IB. 21) arises from the part of parameter space where Z)j_ij-+i, the coefficient of 
Pj_ij-i_i, is small. This region may be identified by parameterizing [66j 

l3i = PiCi, = PiCi, 7i = P2C2, 72 = P2C2 (B.5) 
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2^3 



X2 



X. 1 



X4 



3^3 



X2 



X4 



Xi Xi 

(i) (ii) 

Figure 5: (i) Vertical ladder diagram and (ii) vertical ladder with H-shaped insertion (with numerator 
factor indicated by a dotted line) whose leading log contributions are computed in this appendix. 

where Q runs from to 1, and Q ^ 1 — d. In terms of these parameters, the main contribution 
to the subintegral comes from the region where pi and p2 are both small. Retaining only the 
lowest order terms in A and 

A = aoai + aoa2 + aia2 + 0{pi), (B.6) 
D24 = aoaia2, Dj^ij+i = pip2{ai + 02(1(2 + aoCiG) + 0{p^), 
Dj-1,2 = piao{ai + a2(i) + 0{pl), = P2ao(«i + "2(2) + C'(P2)> 

Dj-1,4 = pia2{ai + ao(i) + 0{pl), -Dj+1,4 = P2a2(ai + ^oG) + C'(P2) 
we may approximate the term in the denominator of eq. (1B.2P as 

D + m^aA ^ Topips + Tipi + r2P2 + Tg (B.7) 

with 

To = («! + ^2(1(2 + aoCiC2)^'j-ij+i 
Ti = ao(ai + a2(i)Pj-i,2 + a2(«i + «oCi)-Pj-i,4 
= ao(ai + a2C2)-Pi+i,2 + "2(01 + aoC2)-P;+i,4 
Ts = Q;oai«2^ + («o«i + c>^o«2 + ttitt2)(«o + ct2)^ • (B-8) 

In the large s hmit, Fq ^ Fi ~ r2 ~ F3. Inserting eq. (\B.7\i into eq. flB.2p . we may approximate 
the integral as 



'double box 



/ dttQ daida2 6 {uq + ai + a2 — 1) {aoai + aoa2 + 0:10:2) / d(id(2J (B. 
Jo Jo 



9) 
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with 



J= dp, dp,- ^^"^"V ^^.3 (B.IO) 

where the factors of pi in the numerator result from the Jacobian in the change of variables. To 
capture the dominant behavior of the integral, we only need integrate pi and p, over a small 
region near the origin; the exact values of the upper limits 771, 772 are unimportant but are both 
taken to be -C 1 to justify the neglect of higher order terms in pi in the integrand and the 
dropping of 5i from the argument of the Dirac delta function. The integral over pi yields 

J^-p-^ — ^log ^ — , ^ / - T^oT^ log To (B.ll) 



(ror3-rir2)2 "V (rir/i + r3)(r2r/2 + r3) ; vlv^ 



where in the last step we used the approximation Fq ^ Fi ~ r2 ~ F3. Inserting eq. (IB. lip into 
eq. (IB.9p . we find 

(B.12) 



^ K\v) logP,_i,,+i 

'double box ~ , ^2 



where v = m? /t and 



K (v) = / aanda^ da-?. / dudG- — - . ^ , : 



daodaida, / d(id(2-, . . ^ > no r / n/ n n 

/o Jo (ai + a2CiC2 + «oCiC2j Foaia2 + (aotti + aoa2 + ai«2j(ao + tt2j'yj 

(B.13) 

Thus in the s — )■ oo limit, integrating out the double box is equivalent to inserting eq. flB.121) 
into the remaining integral. 

We observe that K'{v) may be identified as the coefficient of (logs)/(s^t) in the two-loop 
horizontal ladder diagram in the asymptotic s — ?• oo regime. When v is small, K'{v) can be 
explicitly evaluated to give [30] 



K'{v) = --log^v --TT^hgv + 0{v) . (B.14) 
3 3 

Next we consider an (L — l)-loop vertical ladder integral lL-i,a, supplemented with a numer- 
ator factor of xisX24^^ to make it dual conformal invariant. Consider just the subintegral over 
the jth loop 

T A ^ ^ (5(ao + ai + +71 - 1) .-o .^ 

/single box =/ daodaidPid-ji—— — -— (B.15) 

Jo (roPi7i + Ti/?! + F271 + F3)^ 

where 

Fo = -Pj-ij+i 

Fi = ao -Pj-1,2 + «i Pj-iA 

F2 = Qq -Pj+1,2 + Oil -Pj+1,4 

F3 = aoOit + m^(ao + «i) • (B.16) 

As before, the leading s — )■ 00 behavior of the vertical ladder integral comes from the region of loop 
momentum space where Pj-ij+i is much larger than the other Pik, i.e., where Fq ^ Fi ~ F2 ~ F3. 
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This in turn implies that the main contribution to the subintegral (IB.lSp arises from the region 
of parameter space where /3i and 71 are both smaU. We may therefore approximate eq. (IB.lSp as 



-'single box '~ 

The integral over /3i and 71 yields 



dao dai6{ao + ai — 1) I d(5i I d^i — -— ttt? • (-^■1'^) 







(ro/3i7i + ri/3i + r27i + r3)2 



Tors 



-TiFa V {TiT]i + T3){T2r]2 + T3} J T0T3 



where in the last step we have used the approximation Fq ^ Fi, F2, F3. Hence we obtain for the 
subintegral 

/single box-— „ where K{v) = / daoda^ ■ . (B.19) 

t Pj-i,j+i Jo aoai + V 

Thus, in the s — > 00 limit, integrating out a single box is equivalent to the insertion of eq. fIB.lQP 
in the remaining integral. Here K{v) is just the coefficient of (logs)/(st) in the one-loop box 
diagram in the asymptotic s — > cxd regime. When v is small, it can be explicitly evaluated to give 

K{v) = -2logv + 0{v). (B.20) 

We note that eqs. flB.12p and fIB.lQP are quite similar, except that the powers of Pj_i j+i in 
the denominators differ. This reflects the fact that the two d-lines (described in sec. 12. 2p for 
the double box have length two whereas those for the single box only have unit length. Recall, 
however, that dual conformal invariance implies that the Ilh diagram comes with a factor of 
Pj_ij+i in the numerator (IB. 10 . which cancels one of the factors in the denominator of eq. f lB.12p . 
"promoting" the d-lines to length one. (This in turn raises the number of powers of logs in the 
asymptotic behavior of the integral by two.) 

Taking into account all the numerator factors, we see that the quotient of Ilh and Ii-i^a is 
given by 

= (B.21) 

/i-i.o K(v) ^ ' 

We know however from ref. [301I5T] that 



lim h-i,a = Tr^,K{vf-^ log^-i u + --- (B.22) 



from which we conclude that 



lim hu = \r^^\l K{vf-^ K\v) log^'^ u + --- (B.23) 



which is used in the main body of the paper in eq. ( 12.260 . 
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C Five-point integrals 

In this appendix, we give small-m^ expansions for the Higgs-regulated one- and two-loop five- 
point integrals shown in fig. |3l 

The one-loop box integral I^^^ has the small expansion 



sis^I^^'' = log^ — 2 log ( — ^ ) log 

V S3 / 



m 



+ 21ogSilogS5-log2s3-2Li2(^l-^j _2Li2(^l-^j -y. (C.l) 

Assembling everything into dimensionless ratios Si/rn?, the result can be stated more succinctly 
as 



,,srjr « 21og ^^J log ^-L) - log^ l^-L)-2U,[l-fJ-2U..^l-^J-- + Oirn'). 

(C.2) 

The double-box integral if^"" has the small m? expansion 
s.sll^^)'^ = - ^ log^m^ + log (^^^ log^ 



+ 
+ 



1 2 2 

log m 



3 / S1S2 \ TT 

log 



-^log^ +41og £1 LiJl-^ 



+41og Li2 1 - ^ - 4Li3 1 - ^ + 8i7on 1 - ^ ) - 4(3 

+ C(log° m^), 



logm^ 
(C.3) 



where we use the harmonic polylogarithm function j67] 



iJoii(l -x) = -log(l -a;) log^a; loga;Li2(x) - Li3(x) Cs- 



(C.4) 
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Finally, the pentagon-box integral Jg^''^ has the small expansion 



log(SiS2S3S4) - - log(s2S3S5) 



log^ 



+ 



log Si - 2 log S2 log Si - log S3 log Si + 2 log S4 log Si 



+ log S5 log + ^ log^ S2 + log^ S3 - log^ + ^ log^ S5 
+ log S2 log S3 - 2 log S2 log S4 - 2 log S3 log S4 

-2 log S2 log S5 - 2 log S3 log S5 + log S4 log S5 



Si 



S3 



-Li2 1-- + Li2 1-- -2Li 



S2 



S4 



S2 



oil 



£3 

Si 



1 - 



— 8-ffoii 
+4Li3 ( 1- 

+2 log 



S2 
S5 
S2 

( S1S4S5 \ 



+4 log 



S1S5 



-4Li3 
+ 2 log 
Li2 



1 - 



1 - 



S4 



S5 

- 6iJoii 



- 1 -2Li2 ( 1-- 

S4 



s-i 



£3 
■51. 

/ SlS2S3\ 



■51, 

+ 4Li3 ( 1 



1 - 



S2 



S5 

- 4//0II 



2 ^2 



log m 



£3 
S5 



£4 

S2 



-4Li, 1 - 



£3 
S5 



S2 



Lio 1 - 



£1 

S3 



+ 4 log 



Lio 1 



£1 

S4 



S3 



1 - — -41og 

S4 

S3 



\s4si 



Li2(l-^ 

S5 



Li2 1 + - log si + 2 log S2 log si - 6 log S4 log si 

S5 



- log S5 log^ Si + log^ S2 log Si + 5 log^ S4 log Si - log^ S5 log Si 

1 17 

-2 log S2 log S4 log Si + 2 log S3 log S5 log Si - -TT^ log - Y 

- log^ S3 + 3 log^ S5 - log S2 log^ S3 - 4 log S2 log^ S4 + log S3 log^ S4 

2 

-8 log S2 log^ S5 - 2 log S3 log^ S5 - log S4 log^ S5 - -TT^ log S2 

2 

-2 log^ S2 log S3 - -TT^ log S3 + 6 log^ S2 log S4 + 2 log S2 log S3 log S4 



1 



-TT^ log S4 + 9 log^ S2 log S5 + 2 log^ S3 log S5 + 2 log S2 log S3 log S5 



+ -7r' logS5 + 12C3 



logm^ + 0(\og^ m^^ 



(C.5) 



When we sum over cyclic permutations to obtain the full one- and two-loop amplitudes, the 
polylogarithm functions cancel, resulting in the relatively simple expressions (I3.3p and (I3.8p . 
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